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Abstract 

We study a least squares estimator 9t for the Ornstein-Uhlenbeck 
process, dX t = 9X t dt + adBf , driven by fractional Brownian motion B H 
with Hurst parameter H > i . We prove the strong consistence of 9t (the 
almost surely convergence of 9t to the true parameter 9). We also obtain 
the rate of this convergence when 1/2 < H < 3/4, applying a central limit 
theorem for multiple Wiener integrals. This least squares estimator can 
be used to study other more simulation friendly estimators such as the 
estimator 9t defined by (14. 



1 Introduction 

The Ornstein-Uhlenbeck process X t driven by a certain type of noise Z t is 
described by the Langevin equation 

X t =X -9 [ X s ds + aZ t . 
Jo 

If the parameter 8 is unknown and if the process (X t , < t < T) can be observed 
continuously, then an important problem is to estimate the parameter 9 based 
on the (single path) observation (X t , < t < T). When Z t is the standard 
Brownian motion, this problem has been extensively studied (see for example 
01 [10] and the references therein). The most popular approaches are either 
the maximum likelihood estimators or the least squares estimators, and in this 
case they coincide. Other type of noise processes have also been studied. For 
example, when Z t is an a-stable process maximum likelihood estimators do not 
exist and other approaches are proposed in [5] and [6]. 
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In this paper we study the parameter estimation problem for the Ornstein- 
Uhlenbeck process driven by fractional Brownian motion with Hurst parameter 
H 



X t = X -6 X s ds + aBF, (1.1) 



where 9 > is an unknown parameter. Although the Ornstein-Uhlenbeck pro- 
cess is defined for all H E (0, 1), we assume H > 1/2 in this paper. In [5], the 
the maximum likelihood estimator 9t for the parameter 9 is obtained and has 
the following expression 



9 T = -\J Q 2 {s)dw%\ J Q{s)dZ, 



where 



h,(t,8) = K£ a *- H {t-a)*- H , KH = 2HT[~-H)r[H + ~ 
x - h 2-2H. > 2HT(3-2H)r(H+±) 

<J t — Att t , Aff = 



d r 

Q(t) = — w / k H {t,s)X s ds, 0<t<T: 
dw t J 



Z t = ! k H (t,s)dX s 
Jo 



It is proved that limT^oo 9t = 9 almost surely. 

In this paper we propose two different estimators for the parameter 9 and we 
study their asymptotic behavior. First we introduce an estimator of the form 



e T = e- J^ tdB " , (1.2) 

foXfdt 

where J Q T X t dB^ is a divergence- type integral (see [T], [3J, @], [7] and the refer- 
ences therein), and we call it the least squares estimator. This is motivated by 
the following heuristic argument. The least square estimator aims to minimize 

\X t + 9X t \ 2 dt, 



and this leads to the solution 

a ft**™, 



T 

\ , I I X , 

t „„ .. ■ (1-3) 



So x ? dt 

If H = i, then the integral X t dX t is an Ito stochastic integral which can be 
approximated by forward Riemann sums. However, for H > i the numerical 
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simulation of the estimator 9t seems extremely difficult. For this reason, in this 
case we introduce and study a second estimator 9t, defined in (|4.ip . 

We prove the almost sure convergence of the estimator 9t to 9, as T tends 
to infinity, and derive the rate of convergence, obtaining a central limit theorem 
in the case H S [5,!)- ^he P ro °f °f the central limit theorem is based on 
the characterization of the convergence in law for multiple stochastic integrals 
using the techniques of Malliavin calculus, established recently by Nualart and 
Ortiz-Latorre in [T^] . Finally, we derive the rate of convergence of the estimator 
9t from the rate of convergence of 9t- 



2 Preliminaries 

In this section we first introduce some basic facts on the Malliavin calculus for 
the fractional Brownian motion and recall the main result in [12] concerning the 
central limit theorem for multiple stochastic integrals. 

The fractional Brownian motion with Hurst parameter H £ (0, 1), (Bf 1 , t E 
R) is a zero mean Gaussian process with covariance 

E(5f Bf ) = R H (s, t) = ~ (|t| 2ff + \s\ 2H \t s\ 2H ) . (2.1) 

We assume that B H is defined on a complete probability space (£l,A, P) such 
that A is generated by B H . Fix a time interval [0, T]. Denote by S the set of 
real valued step functions on [0, T] and let 7i be the Hilbert space defined as 
the closure of £ with respect to the scalar product 

(l[o,t]> l[o,s])w = Rn(t, s), 

where Rr is the covariance function of the fBm, given in (|2.1[) . The mapping 
l[o,t]' — * can be extended to a linear isometry between Tt and the Gaussian 
space Hi spanned by B H . We denote this isometry by ip 1 — ► B H (ip). For 
H = \ we have H = L 2 ([0, T}), whereas for H > \ we have L^([Q, T]) C H and 
for ip, ip £ Ltt([0, T}) we have 

{ip,iP) n = a H [ [ Vsipt\t-s\ 2H - 2 dsdt, (2.2) 
Jo Jo 

where a H = H(2H - 1). 

Let iS be the space of smooth and cylindrical random variables of the form 

F = f(B H (^),...,B H ^ n )), (2.3) 

where / S C£°(IR n ) (/ and all its partial derivatives are bounded). For a random 
variable F of the form (|2.3|) we define its Malliavin derivative as the W-valued 
random variable 

^ = E|^(^i),---,^ n ))^. 
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By iteration, one can define the mth derivative D m F, which is an element of 
L 2 (tt;H® m ), for every m > 2. For m > 1, B" 1 ' 2 denotes the closure of S with 
respect to the norm || • || m ,2, defined by the relation 

m 

\\F\\ 2 m ,2 = E[|F| 2 ]+^E(||^||^)- 

i=l 

Let 5 be the adjoint of the operator D, also called the divergence operator. A 
random element u G L 2 (il,H.) belongs to the domain of <5, denoted Dom(<5), if 
and only if it verifies 

\E(DF,u) H \<c u \\F\\ L 2, 

for any F £ D ' , where c u is a constant depending only on u. If u € Dom(<5), 
then the random variable 5(u) is defined by the duality relationship 

E(F8(u)) = E(DF,u) H , (2.4) 

which holds for every F <G D 1 ' 2 . The divergence operator 5 is also called the 
Skorohod integral because in the case of the Brownian motion it coincides with 
the anticipating stochastic integral introduced by Skorohod in [T5|. We will 
make use of the notation S(u) = u t dB^ . 

For every n > 1, let TC n be the nth Wiener chaos of B, that is, the closed lin- 
ear subspace of L 2 (fi, A, P) generated by the random variables {H n (B H (h)) , h S 
H, \\h\\-n — 1}, where H n is the nth Hermite polynomial. The mapping I n (h® n ) = 
n\H n (B H (h)j provides a linear isometry between the symmetric tensor product 
7i 0n and Ti n . For H — \,I n coincides with the multiple Ito stochastic integral. 
On the other hand, l n (h® n ) coincides with the iterated divergence 8 n {h® n ). 

We will make use of the following central limit theorem for multiple stochas- 
tic integrals (see [12]). 

Theorem 2.1 Let {F n , n > 1} be a sequence of random variables in the p- 
th Wiener chaos, p > 2, such that limjj^oo E(i 7 ' 2 ) = a 2 . Then the following 
conditions are equivalent: 

(i) F n converges in law to N(0,a 2 ) as n tends to infinity. 

(ii) WDFnW 2 , converges in L 2 to a constant as n tends to infinity. 

Remark. In [12] it is proved that (i) is equivalent to the fact that H-D-Fnll^ 
converges in L 2 to pa 2 as n tends to infinity. If we assume (ii), the limit of 
H-DFnllf^ must be equal to pa 2 because 

E(|| J D^|| 2 ,)=pE(^ 2 ). 
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3 Asymptotic behavior of the least square esti- 
mator 

Consider Equation driven by a fractional Brownian motion B H with Hurst 
parameter H > |. Suppose that Xq = and 9 > 0. The solution is given by 



X t = a /V 9 «->dB?, 
J o 



(3.1) 



where the stochastic integral is an Ito integral if H = h and a path-wise 

Riemann-Stieltjes integral if H > i. Let 9t be the least squares estimator 
defined in (|1.2|) . The next lemma provides a useful alternative expression for 

Lemma 3.1 Suppose that H > |. T/ien 

fe— + (3.2) 

Proof Using the relation between the divergence integral and the path-wise 
Riemann-Stieltjes integral (see Theorem 3.12 and Equation (3.6) of [3]) we can 
write 

fT r T rT ft 

salt 



f X t odB? = f X t dB? + a H [ f D s X t {t~ s) 2H - 2 ols 
Jo Jo Jo Jo 

= [ T X t dB? + aa H [ T f e~ 8 ^ (t - s) 2H - 2 dsdt 
Jo Jo Jo 

= f XtdBf + aa H [ [ £ 2H ~ 2 e- e td£dt . 
Jo Jo Jo 



As a consequence, we obtain 

rT 



e T = e-a^-^ t - + a 2 Jo Jo z, — . 3.3 

fixfdt tfxfdt 

On the other hand, 

T T T T 

a f X t ° dB? = I X t o dX t + 9 [ Xfdt = -X% + 9 1 X 2 dt. (3.4) 
Jo Jo Jo 2 Jo 

Substituting §3$ into Q yields (021). □ 

The next theorem establishes the strong consistency of this estimator. 

Theorem 3.2 If H >\, then 

9 T ^9 (3.5) 
almost surely, as T tends to infinity. 
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In order to prove this theorem we make use of the following technical result. 
Lemma 3.3 Assume H > |. Then, 

1 r T 

- / X 2 dt -> a 2 0- 2H HT(2H), (3.6) 

2 Jo 

almost surely and in L 2 , as T tends to infinity. 
Proof For every t > define 

r t = a f e' e ^- s UB^ =X t + e~ 6 % (3.7) 

J — oo 

where £ = a e 9s dB^ . The stochastic process (5^,t > 0) is Gaussian, sta- 
tionary and ergodic. For H = ^ this is well-known and for H > ^ this is proved 
in [2J. Then, the ergodic theorem implies that 



f 



T 

2 



as T tends to infinity, almost surely and in L . This implies that 

i-T 

X 2 dt -> E(Y 2 ), 



1 

as T tends to infinity, almost surely and in L 2 . If H = i, we know that 
E(y 2 ) = f^, which implies ([Bl]) , If # > i, using yields 

p OO /* OO 

E(X 2 ) = a H a 2 / / e~ 9 ^ \u - s | 2ff - 2 du<fa, 
Jo Jo 

and (|3.6p follows from Lemma ISTTl □ 

Proof of Theorem \3.2\ In the case ff = |, taking into account that the process 

X s dB s ,t > 0^ is a martingale with quadratic variation J* X 2 ds it follows 

that 9t — > ^ almost surely, as T tends to infinity. 

Now let -ff > 1/2. From Lemma 15721 we deduce that almost surely 

km = o. (3.8) 

It is easy to check that this convergence also holds in L 2 . Then we conclude the 
proof using Lemma [331 <|3.8j) . and 

lim i / T f e H ^e-°^dt = B l -^T{2H-l). 

□ 

The next theorem provides the convergence in distribution to a Gaussian 
law of the fluctuations in the almost sure convergence (|3.5j) . 
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Theorem 3.4 Suppose H € [5, |). Let (X t ,t E [0,T]) be given by \3.1\l , then 

VT \e T - e\ 4. N(0, 6a 2 H ) , (3.9) 
as T tends to infinity, where 

Proof We have 



j T x 2 dt ffxfdt f'Xfdt 



(3.11) 



F T = ^/ 2 (e- e l*- s l). (3.12) 



where Ft is the double stochastic integral 

r 2 

2Vr 

From Lemma T3. 31 wc know that A Xfdt converges almost surely and in L 2 , 
as T tends to infinity to a 2 9^ 2H HT{2H). Then, it suffices to show that F T 
converges in law as T tends to infinity to a centered normal distribution. In 
order to show this convergence we will apply Theorem l2.1l to a given sequence of 
random variables in the second chaos Fr k , where | 00 as k tends to infinity. 
To simplify we assume that T — 1,2,.... The proof then follows from the 
following facts: 

(i) E(i*y) converges to (9 1_4fl a A o~H, where 

fa = H^H - WW + 

as T tends to infinity. 

(ii) ||Z)Fy|||j converges in L 2 to a constant as T tends to infinity. 

Step 1: Proof of (i) Suppose first that H = In this case, by the isometry 
of the Ito integral we obtain 

E(F 2 ) = C r f e~ 26{t - s] dsdt ' T ' ' T 



T J J T \ 26 46> 2 



which implies that 



a 4 



A,*™ 29 



This implies the desired result because Si 
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Now, let H £ (|, |). In this case, by the isometry property of the double 
stochastic integral I2 , the variance of Ft is given by 

EM = ^It, (3-13) 



where 



I T = I e- e ^^- e ^^\8 2 -8^ 2H - 2 \u2-ux\ 2H - 2 dud8. (3.14) 



'[0,T] 4 

By Lemma 15.31 in the Appendix we obtain that 

lim EOF!) - e^^a^H- 

T — >oo 

Step 2: Proof of (ii) For s < T we have 

D s F T = ^ + 4= I' 'e-W-)dB!' . 



Suppose first that H = \. In this case, 



2\ 



\\DF T \\ 2 n = ^ £^X 2 + 2aX s j\~ e ^dB t + a 2 {j~ e-°^dB t \ \ ds 

- 4« 4. a {2) 4- /l (3) 

We already know from (|3.6[) that converges in L 2 to ^ as T tends to 
infinity. The third term can be written as 

^ - £ jf (f-— - £ f (i*- s< - ■>«.)'*- 

so it also converges in L 2 to |s a T tends to infinity. Finally we can show that 

lim E((4 2) ) 2 ) = 0. 

In fact, we have 

E((4 2) ) 2 ) = % J { <t} E (x s X u ^j\- 6 ^dB t j U\- e ^dBMdsdu 

= ^ I ( f S e-^ s+u - 2 ^dA ( C ' e-^—^dt] dsdu 

8ct 8 



{s<u<T} \J0 



4^ 2 T 2 J {s<u < T} 



(e 2f>s -\){e- 2eu -e- 2eT )dsdu, 
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which clearly converges to zero as T tends to infinity. Therefore, 
converges to ^- in L 2 . 

Suppose now that H > |. From (12. 2ft we have 

WDF T f H = [ T I" ( Xs + a f T e-^dB? 



o Jo 



\u — s 



2H - 2 duds 



We have to prove that ||Z?Ft||^ converges to a constant in L 2 as T tends to 
infinity. In fact, 



2 i-T ,-T 

\DF T \\ 2 H = 



^ jf (x s X u + 2aX u jf e- 0(t - s UB? 



-a 1 / , " ' -,//>'/' / e-W-^dBf ] ( /-.s-| 2 " 2 ,/ (( </.s 

. ~(2) . ~(3) 



T 



For the term , since X t is Gaussian we can write 



E (|4 1} - H^t^I 2 ) = 2 / E (X s X t ) E (X U X„) 

>■ / /[n T14 



'[0,T] 4 

x |u - s| 2H " 2 |w - t\ 2H ~ 2 dudvdsdt. 



By Lemma [5741 

- / E (XrX t ) E (X U X H )(T - u) 2H ~ 2 \v ~ t\ 2H - 2 dudvdt 
< I / ( T - t) 2H - 2 | M - V | 2ff - 2 (T - U ) 2ff - 2 |t; - t\ 2H - 2 dudvdt 

* J[0,T\* 

<C 2 gH T 8H - 6 [ (l-t) 2H - 2 \u-v\ 2H - 2 (l-u) 2H - 2 \v-t\ 2H - 2 dudvdt, 
J[o,i] 3 

which converges to as T tends to infinity when H < |. Hence, by l'Hopital 
rule, E — E(Cy 1 ' ) )| 2 ^ converges to as T tends to infinity. In the same 

way we can prove that E ^|C^ — E(C^'')| 2 ^ converges to zero as T tends to 
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infinity, for i = 2, 3 when H < 3/4. By triangular inequality, we see that 
E[(\\DF T \\ 2 H -E(\\DF T f H )y 

= e (\cp + 4 2) + 4 3) - e(4 1} + 4 2) + 4 3) )i 2 ) 



< 



9^£;(|4 l) -E(4 4) )| 2 



0. 



Taking into account that 

Km E(||I>F T ||*) = 2 Km E(F|), 

1 — >oo 1 — >oo 

we conclude the proof of (ii). This completes the proof of the theorem. □ 

If one replaces the ltd type integral in (|1.2p by the path-wise Riemann- 
Stieltjes integral, then we can obtain the following estimator 



/o X t o dX t 



f'Xfdt 2fJX?dt 
which converges to zero in L 2 as T tends to infinity from Lemma l3.3l and (|3.8p . 



4 An alternative estimator 

Suppose in this section that H > |. We introduce the following estimator 

i 

/ 1 



a 2 HT(2H)T 



X 2 dt 



(4.1) 



From p.6p . we see that 9t converges to 8 almost surely as T — > oo. Theorem 
allows us to derive the rate of convergence in the approximation of 6 by 9t- 



Theorem 4.1 Suppose H £ (±, |). Then 



Vt 



9 t -9 



N 



(2H) 



2 

2 a H 



as T tends to infinity, where o~h is defined in h3.10\ ). 
Proof From Equation (|3.2j) . we have 

' X 2 dt = <p«h£ !te H -y-*di-dt-x$i2 



(4.2) 



10 



Thus 



- VT 



HT(2H) 



a H ±£tie"-*e-<Xdt;dt-^ 



q 2H 

7 T 



From Lemma [5721 it follows that 



T r t 



JO 



2H-2-6 



X, 



d£dt - ^ = <*hT(2H - l)6 x - 2H + o(^= 



T 



where o{-^=) denotes a random variable Ht such that VtHt converges to zero 
almost surely as T tends to infinity. Therefore, 



HT(2H) 



aH±fotiZ 2H - 2 e-<>ZdZdt-^ 



On the other hand, we can write 



(4.3) 



Q2H _ Q~2H 



= v T 



2 [ ° T 



where 9^ is a random point between 9 and Ot- From Theorem 13.41 we obtain 
the following convergence in law as T tends to infinity: 



2H 



N 0, 



1 



{2Hf 



> 2H a 2 H 



(4.4) 



Finally, from the decomposition 

Vt 



HT(2H) 



q 2H 
7 T 



and using (|4. 3[) and (|4.4p we deduce the desired convergence. 



□ 



5 Appendix 

In the sequel we present some calculations used in the paper. 
Lemma 5.1 For any H £ (i l) 



/■OO /-00 

(2H-1) / 
Jo Jo 
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Proof We can write, by the change-of-variables u — s = 



OO pU 



o Jo Jo Jo 

poo pit 

= 2 / e- 2u+x x 2H - 2 dxdu. 
Jo Jo 

Integrating first in the variable u and using that (2H — 1)T(2H — 1) = T(2H) 
we conclude the proof. □ 

Lemma 5.2 Let Y t be the stationary Gaussian process defined in \3. 7| ), where 
H > 1/2. Then, for any a > 0, converges almost surely to zero as T tends 
to infinity. 



Proof The covariance of the process is Yt is, using Lemma 15.11 to compute 
Var(0, 



Cov(Yo,r t ) = e~ 6t ^U Z + a J e 6u dB* 

Vw(0+o- 2 a H f [ 
Jo J-; 



9u+6v 



\u-v\ 2H ~ 2 dudv 



= e 



-01 



-01 



i>t poo 

Var(0 + o*a H / / e- 9x+26v x 2H - 2 dxdv 

Jo Jv 



a 2 6- 2H HT(2H) + a 2 { 9 1 - 2H HT(2H)t - -t 2H + o(t 2H ) } 



a 2 e~ 2H HT{2H) 



q2H 



a 2 T{2H 



T f H +o{i 2 ")} 



Then the result lemma from Theorem 3.1 of Pickands [T^] . 
Lemma 5.3 Let Lt given by pJ^ ). When \<H<\, we have 



iim * = e 1 -^^, 



wher 



□ 



(5.1) 



lH = (8H - 2) T{2H - l) 2 + (16-ff - 4) 



T(2H - l)r(3 - AH)T{AH - 2) 



r(2-2F) 

Proof Taking the derivative with respect to T we have 



dI T 

It 



= 4 / e- e ^ T - u ^- 6 ^-^ {T - Sl ) 2H - 2 \u 2 -u l \ 2H - 2 du 1 du 2 ds 1 . (5.2) 

J{0,T] s 
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dI T 
It 



Making the change of variable T — u 2 = x\ , T — u% = x%, and T — s% = X3 yields 
4 / e-ozx-ete-* sl^-2| Xi _ X2 fH-2 dxidx2dx ^ 

J[0,T] 3 

As a consequence, 

rim ^ = 4 / e -«x 1 -©|x B -« 8 | a; a^-2| a . i _ ^-2^^^ 

and this integral is finite. Indeed, we can decompose this integral into the 
integrals in the six disjoint regions {x a m < ^ CT (2) < x <?(3)}> where a runs over 
all permutations of the indices {1,2,3}. In the case x\ < x 3 < x 2 making the 
change of variables X\ — a, x 3 — X\ — b, and x 2 — x 3 = c, we obtain 

e -0(a+c )(a + b) 2H-2 {b + c fH-2 dadMc 

[0,oo) 3 

< / e'^ a+ % AH - A dadbd Cl 

^[0,00)3 

which is finite because H < §. The other cases are simpler and can be handled 
in a similar way. We can write 

f e -e Xl -e\ X2 - X3 \ x 2H-2\ Xi _ X2 fH-2 dj , idx2dx3 = e l-4H d[{ ^ 

J[0,Oo) 3 

where 

d H = f e- x - ly - zl z 2H - 2 \x-y\ 2H - 2 dxdydz. (5.3) 

J[0,oo) 3 

The integral in (|5.3p can be simplified as follows. First we make the change of 
variables y 1— > w, where w = y — x, and we obtain 

/>oo />oo />oo 

d H = / / e- x -\* +w - z \z 2H - 2 \w\ 2H - 2 dwdxdz 

JO J0 J-x 

= / / / e -( 2 -+— z )z 2H - 2 \w\ 2H - 2 dwdxdz 

Jo Jo J z-x 

pOO P OO pZ — X 

e-^z 2H - 2 \w\ 2H ~ 2 dwdxdz. 



Jo 



Integrating in x we get 

1 />oo />oo 

2 Jo »/-oo 

[{z -w)- ((-to) V 0)], e- (z - w V H - 2 |™| 2i? - 2 dujcfe. 



OO /"OO 
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Therefore, 



1 

2 

+: 
+ 



OO pOO 



IQ JO 

poo poo 



JO 

oo poo 



JO 

oo poo 



JO 



e- 2 ^- w ^°^ w - z h 2H - 2 w 2H - 2 dwdz 

e -(*+™) z 2H-2 w 2H-2 dwdz 

[(z - w)} + e-( z -^z 2H ~ 2 w 2H ~ 2 dwdz 
-^z^-'w^^dwdz, 



and we obtain 



where 



d H = f H +[2H-^)T(2H-l) 2 , 



(5.4) 



OO fZ 



f H = (l + z- w)e-( z - w h 2H - 2 w 2H ~ 2 dwdz. 

Jo Jo 

Making the change-of-variables z ~ w = x yields 



oo poo 



fa 



(1 + x)e-*(w + xY H - z w 2H - I dwdx. 



JO 



Substituting the equality {w + x) 2R - 2 = T(2 \ H) f °° e~^ w+x ^ 1 - 2H d^ in f H 
we obtain 



-I P OC 1 p OO P OG 

f H = — n om / / / (l + x)e-*-e^v? B -*e- 2H dtdwdx 
- Zti) Jo Jo Jo 

{l + x)e- x -^ x ( 1 - AH didx 

y 

(1 + x)e- x x iH - 3 dx 



r(2- Jo j j 
T(2H-l) pCO r °° 



T(2-2H)J Jo 
_ T(2H- l)r(3-4ff) 
_ r(2-2ff) jo 

iN r(2g-i)r(3-4g)r(4g-2) 
= (4ff ~ 1} r(2-2ff) • 

Finally from (|5.4[) and (|5.5p we get the desired result. 
Lemma 5.4 Let X t be given by U.2\) . We have 



E 



e- 6 ^- s UBf / e - e ^dB H 



< C e , H \t-s 



2H-2 



and 



E[X t X s ]<a 2 CgM\t-s\ 2H - 2 , 
for some constant Ce,H > 0. 



(5.5) 
□ 

(5.6) 
(5.7) 
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Proof Let us assume that s < t. We can write using 



E 



T i-T 



1{t-s) p -8{ri-t)\ 



\C-v\ 2H - 2 d^dr 1 = a H (B^+B^) 



where 



and 



R (i) _ 



T i-T 



Ki-s) e -e( v -t) 



o(2) _ 



t Jt 

T r t 



t J s 



e " 'e 
e -9tt-) e -fl(-J-*)| € _^|aff-a^. 



It is easy to see that is bounded by C^jje e '* s '. The second term can be 
estimated as follows 



B, 



(2) 



/* ( T ? e- e tt- s+ y + t-Vy 2H - 2 dydt 

Js Jt-f 



T-s 







(T-w)At 



(t-y)Vs 



20 



20 



T-s 



< _ / e -8(y+2(t-y)Vs-s-t) y 2H-2 dy 



T-s 



-8(y+s-t)„2H-2 



1 



2(9 



e -e(y + s-t) y 2H-2 dy 



< C e , H \t~s 



2H-2 



2H-2 



T-s 



e -e(y+s-t) dy + C e / y 2H - 2 dy 

t-s JO 



This proves (|5.6|) . The inequality (|5.7|) can be proved in a similar way (see also 
0). □ 



References 



[1] Biagini, F.; Hu, Y.; 0ksendal, B. and Zhang, T. Stochastic calculus for 
fractional Brownian motion and applications. Springer, 2008. 

[2] Cheridito, P.; Kawaguchi H. and Maejima, M. Fractional Ornstcin- 
Uhlenbeck processes. Electronic Journal of Probability 8 (2003) 1-14. 

[3] Duncan, T. E.; Hu, Y. and Pasik-Duncan, B. Stochastic calculus for frac- 
tional Brownian motion. I. Theory. SIAM J. Control Optim. 38 (2000) 
582-612. 



15 



[4] Hu, Y. Integral transformations and anticipative calculus for fractional 
Brownian motions. Mem. Amer. Math. Soc. 175 (2005) 825. 

[5] Hu, Y. and Long, H. Parameter estimation for Ornstein-Uhlenbeck pro- 
cesses driven by a-stable Levy motions. Commun. Stock. Anal. 1 (2007) 
175-192. 

[6] Hu, Y. and Long, H. Least Squares Estimator for Ornstein-Uhlenbeck Pro- 
cesses Driven by alpha-Stable Motions. Stochastic Process. Appl. In press. 

[7] Hu, Y. and 0ksendal, B. Fractional white noise calculus and applications to 
finance. Infin. Dimens. Anal. Quantum Probab. Relat. Top. 6 (2003) 1-32. 

[8] Kleptsyna, M. L. and Le Breton, A. Statistical analysis of the fractional 
Ornstein-Uhlenbeck type process. Stat. Inference Stock. Process. 5 (2002) 
229-248. 

[9] Kutoyants, Yu. A. Statistical Inference for Ergodic Diffusion Processes. 
Springer, 2004. 

[10] Liptscr, R.S. and Shiryaev, A.N. Statistics of Random Processes: II Appli- 
cations. Second Edition, Applications of Mathematics, Springer, 2001. 

[11] Nualart, D. The Malliavin Calculus and Related Topics. Second edition, 
Springer, 2006. 

[12] Nualart, D. and Ortiz-Latorre, S. Central limit theorems for multiple 
stochastic integrals and Malliavin calculus. Stochastic Process. Appl. 118 
(2008) 614-628. 

[13] Nualart, D. and Peccati, G. Central limit theorems for sequences of multiple 
stochastic integrals. Ann. Probab. 33 (2005) 177-193. 

[14] Pickands, J. Asymptotic properties of the maximum in a stationary Gaus- 
sian process. Trans. Amer. Math. Soc. 145 (1969) 75-86. 

[15] Skorohod, A. V. On a generalization of a stochastic integral. Theory Probab. 
Appl. 20 (1975) 219-233. 



16 



